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1. INTRODUCTION 

With the advent of a new generation of large- 
scale detectors of cosmic radiation, the observa- 
tion of high-energy neutrinos produced in distant 
astrophysical sites (or possibly by other, more 
exotic, mechanisms) has become one of the ma- 
jor challenges of astroparticle physics [1]. Given 
the large distances traveled, approximately equal 
fluxes in , and i^r are expected on Earth as 
a result of flavour mixing and oscillations [2|3| . 

The detection of Earth-skimming Vr's (through 
the observation of the shower induced by the 
emerging r lepton 4 ) has become a very promis- 
ing strategy for the observation of ultra-high en- 
ergy (UHE) cosmic neutrinos. This method has 
been used recently by cosmic-ray experiments 
such as HiRes and the Pierre Auger Observatory 
to put competitive limits on the flux of cosmic 
neutrinos in the energy range 10~^ EeV < < 
102 EeV 

The sensitivity to such UHE Earth-skimming 
neutrinos crucially depends on the parameters of 
the Vt and r propagation through the terrestrial 
crust, and on the correct estimation of the flux 
of r's that emerge from the Earth (see Fig. [1]). 
This problem has been widely discussed in differ- 
ent contexts and with different approximations 
(see e.g. the references provided in |7|8| ). An 
exhaustive treatment should account for r and 
Vt neutral- and charged-current interactions with 



nucleons, r decay and energy losses. 

However the full coupled transport equations 
admit no analytical solution, and they are usu- 
ally treated numerically (or semi-analytically) in 
a simplified framework where several effects are 
neglected, such as the possibility of multiple re- 
generations of the Vr, the weak interactions of 
the T, and the stochastic nature of its energy 
losses. Significant uncertainties in the calcula- 
tion also derive from the poor knowledge of the 
relevant cross-sections at UHE, where no direct 
measurements exist. This is particularly true for 
the neutrino-nucleon interactions and for the pho- 
tonuclear contribution to the tau energy losses, 
which both rely on extrapolations of the structure 
functions beyond the range of (x, Q^) probed at 
accelerators and which could be affected by the 
onset of new physics beyond the Standard Model. 

These considerations demonstrate the need to 
assess the impact of such simplifications on the 
determination of the final flux of emerging r's. 
We present here such a study of the propagation 
in standard rock of tau leptons and neutrinos with 
both mono-energetic and power-law spectra. 

The common framework and the strategy used 
to solve the transport equations is briefly de- 
scribed in Section [21 In Section [3] we focus on the 
r energy losses and compare the results obtained 
with a full stochastic treatment and within the 
continuous approximation. Then, we investigate 
in Section |3] the impact of the ^ t ^ v-^ re- 
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generation chain both in the standard case and in 
non-standard scenarios for the neutrino-nucleon 
interactions and for the tau energy losses. The 
conclusions are presented in Section [51 where we 
also discuss the impact of our studies on the sen- 
sitivity of current extensive air shower detectors 
to Earth-skimming Vr^s. 

2. i^r AND T PROPAGATION: GEN- 
ERAL TRANSPORT EQUATIONS 

The geometry of the propagation problem is 
described in Fig[Tl where an example of regen- 
eration chain through multiple CC interactions 
and T decays is sketched. Given a beam of paral- 
lel neutrinos incident on the Earth, the problem 
becomes uni-dimensional. The flux of the tau lep- 
tons that emerge only depends on the amount of 
crossed rock which, in the hypothesis of a spher- 
ical Earth with a constant density in its crust, is 
given by the incident nadir angle a, which is also 
the angle of the emerging tau. 

Two coupled, integro-diffcrcntial equations de- 
scribe the evolution of the and t fluxes, 
and ^r, along their paths through the Earth, ac- 
counting for all possible production and absorp- 
tion processes taking place within an infinitesimal 
displacement dx: 
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where and X'^^ are the mean free paths 
corresponding respectively to NC and CC inter- 
actions of the incident while a^'~^ corresponds 
to the r CC interaction, which regenerates a v^- 
Ft- is the tau lepton lifetime. Similarly, the equa- 
tion for the T reads: 
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where the index i runs over the r en- 
ergy loss processes, namely pair production, 
bremsstrahlung and photonuclear interactions. 

At the energies of interest for this study, the 
relevant electromagnetic processes that the tau 
lepton undergoes are the bremsstrahlung, the 
pair production and the photonuclear interac- 
tions. The dominant source of t energy losses 
ai Et > Q EeV comes from the photonuclear in- 
teractions, whose description relies on the proper 
modelisation of the nucleon structure functions 
at low X and high . Several calculations of 
the corresponding cross section, based on dif- 
ferent theoretical models, can be found in the 
Iiterature [9ll0lllll2j . All of them give compa- 
rable results except for [llj which gives a signifi- 
cantly higher loss rate, and for which gives a 
rate more than a factor of 2 lower than the stan- 
dard one at the highest energies 10^^ GeV). 

As for the neutrino-nucleon cross-section, a 
benchmark parameterization based on recent 
HERA data is given in [13]. The same paper 
points out another, more speculative, approach 
based on the color glass condensate formalism 
which gives a much lower cross-section. On 
the other hand, plenty of models using new 
physics predict instead an enhancement of <j,yM 
at UHEjig. 




dx 
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Figure 1. Geometry of the transport problem. 
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3. TAU ENERGY LOSSES: CONTINU- 
OUS VS STOCHASTIC APPROACH 

In order to investigate whether the stochastic- 
ity of this process has an effect on the description 
of the propagation of UHE r's through the Earth, 
we will first work in a slightly simplified frame- 
work where we neglect the r weak interactions 
(a reasonable assumption at UHE) as well as the 
possibility of r — > z/,- — > r regeneration chains (an 
effect which will be studied in Section |4]). For the 
T photonuclear interactions, we use here the stan- 
dard parameterization of the PDFs given in (lOj 
based on Regge theory. 

When the differential cross-sections exhibit a 
peak near y=0, as it is indeed the case for the 
processes that govern the r energy losses, the in- 
tegrals are dominated by the behavior of the inte- 
grands around 0, in such a way that an expansion 
of these integrands can be performed. At first 
order in y, and with the simplificating assump- 
tions described hereabove, this gives the following 
transport equation for the t: 

where we have introduced the standard notation 
/3(^) = ^E/^"'""'dyy^(£;,2/). (4) 

Within the approximation of continuous energy 
losses, the average energy lost by the t per unit 
distance di^r/da; is, at UHE, assumed to be pro- 
portional to the mean inelasticity of each process 
which is directly related to (3{E), yielding the fa- 
miliar expression : 

^ = -p P{Er) Er. (5) 

da; 

In that case, Eqn.[3]can be easily integrated, lead- 
ing to the following expression : 

<^{E,x) = $o(^o) X 

xexp£d.(A,(^„)-3-i^) (6) 

where ^{E) = pEf3{E) and E^ the solution of 



In the following, we use Eqn. [6] to compute any 
propagated flux of tau leptons when referring to 
as the continuous energy losses (CEL) approxi- 
mation. 

What we want now is to compare the accu- 
racy of results obtained by using this approxima- 
tion with those from a fully simulated transport 
equation including the stochastic terms. There- 
fore we use a Monte-Carlo generator sampling all 
the interactions, which separates the losses into 
two components jl6| : a continuous one where the 
rate of the losses is large (y G [ymin,ycut]), and 
a stochastic one where the differential cross sec- 
tions lead to more catastrophic losses but with a 
weaker rate {y e [yciLt,ymax])- A good compro- 
mise to reproduce the stochastic features using a 
reasonably fast code is to take j/j. 



- 10^3 



We show in Fig. [2] the simulated distributions 
of the energy of an incident monoenergetic tau 
beam after crossing different depths of standard 
rock, as obtained in the stochastic framework, to- 
gether with the mean of the energy distributions, 
which are the values used within the CEL approx- 
imation. 

After I km, many of the simulated events 
still carry a large fraction of the initial energy 
Eq, meaning that these particles did not un- 
dergo many interactions. However, the distri- 




/os„ (Energy [GeV]) 



Figure 2. Energy logio-distributions of taus oi Eq = 3 ■ 
10^ GeV after propagation in standard rock for 3 different 
depths (1, 5 and 10 km). Also shown are the mean values 
of the energy distributions. 
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bution is asymmetric and there is a large tail of 
events undergoing hard losses. For longer paths 
in rock, fluctuations in the energy losses increase, 
resulting in a broadening of the distribution and 
a smoothening of the high-energy cutoff. The 
stochasticity of r energy losses thus globally leads 
to a smaller survival probability for the r than in 
the CEL approximation; the difference in the r 
range is however always small, as can be seen in 
FigEl 

We have also looked for possible distorsions on 
the propagated fluxes due to the stochastic effects 
in the case of continuous E^^ and E^^ tau injec- 
tion spectra in the range 10* GeV < Er <3- 10" 
GeV. The cutoff for the maximum energy is cho- 
sen to be sharp to exhibit most clearly the differ- 
ent behaviors of the propagated fluxes. We com- 
pare the results obtained with the stochastic and 
the continuous energy losses in Fig [4] for the case 
of an E~^ incident r flux, where the differences 
are more visible. 

Here again, the sharp cutoff present in the CEL 
approximation gets smoothed when accouting for 
stochastic processes: there are indeed fluctua- 
tions affecting a small fraction of particles which 
undergo less interactions and less hard losses. 
The energy range on which this broadening oc- 




Distance \km\ 



Figure 3. Range R{Eo) = dx Psurv{Eo,x) of the 
tau calculated with continuous (top curve) and stochastic 
(bottom curve) energy losses. For comparison, the decay 
length of the tau (dashed line) is also shown. 



curs increases with the depth traversed by the 
tau, and with the hardness of the incident flux. 
It is the only important distorsion induced by the 
stochastic processes. This means that, as far as 
we are dealing with continuous incident spectra 
behaving as power-laws, there is a compensation 
between positive and negative fluctuations of the 
energy losses everywhere in the considered energy 
range, except near the high-energy border. This 
compensation allows to use the continuous energy 
loss approximation as the correct mean value of 
the propagated spectrum. 

4. EFFECT OF i^r REGENERATION 

In this section, we present the results of a study 
of the effects of I'r and r regeneration on the 
flux of emerging r leptons in different scenar- 
ios for the neutrino-nucleon cross-section and for 
the T energy losses. For this, we go back to the 
full propagation equations ([1]) and including 
all possible regeneration terms, except that we 
replace the stochastic energy loss term for the 
T by the CEL expression as given in Eq. ([3]), 
in concordance with the conclusions of Section 
El As stated in Section [21 different calculations 




Energy [GeV] 

Figure 4. Comparison of the emerging t fluxes (normal- 
ized to the injection flux) obtained with the CEL approxi- 
mation (Eqn[6ll and with the full Monte-Carlo calculation 
including the stochastic effects (squared dots). An 
flux between 10* GeV and 3 ■ 10^^ GeV is used at injec- 
tion. Three distances of propagation in rock are shown: 1 
km, 5 km and 10 km. 
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of these quantities exist in the literature. For 
the purpose of our analysis we have extracted 
a panel of ad-hoc parameterizations representing 
respectively a low, a standard (std) and a high 
value of either the t energy loss parameter /3 or 
neutrino-nucleon cross section a^N (the explicit 
expressions are given in [8|). We solve the trans- 
port equations through a Monte Carlo simulation 
(which was also cross-checked against an itera- 
tive semi-analytical method: see the appendix of 
[8|) and compare the results obtained for different 
combinations of j3i ® a^^ {iji—low, std, high). 

As a starting point, one should notice that the 
regeneration is a second order process which re- 
quires the T to be converted back into a Vr before 
loosing too much energy, and then this Ur to un- 
dergo another CC weak interaction to produce a 
T again. Hence, one expects that a higher energy 
loss or a lower neutrino-nucleon cross-section will 
reduce the effect of the regeneration, while, on the 
other hand, a higher uN cross-section or a lower 
energy loss will enhance it. 

To verify this effect, we have first propagated 
through the Earth incident beams of monoener- 
getic J^r's injected at an angle a up to 3°, which 
roughly corresponds to the angular range of sen- 
sitivity of a surface detector such as the one of 
the Pierre Auger Observatory. We find that for a 
standard choice of parameterizations (3std €5 cr^^^ 
(see [8] for more detail), the effect of regenera- 
tion is negligible for E^^ < 10^ GeV, but not at 
higher energies where it can significantly increase 
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Figure 5. Flux of emerging taus corresponding to in- 
cident neutrinos with an angle a up to 3° or 15° and 
^(E^) oc E^^ using the standard combination of param- 
eterisations (Tstd ® Pstd- The regeneration is neglected in 
the dashed histograms and included in the solid ones. 



the flux of emerging r's (15% more r's for 30 EeV 
incident neutrinos). 

As can be seen from FiglSl this effect is washed 
out in the case of an incident neutrino flux oc 
£'~^, for which the contribution of the highest- 
energy part of the neutrino spectrum is negligible; 
but it could affect the spectrum of emerging r's 
if the neutrino flux were harder. Moreover, as 
shown in the same figure, even in the standard 
picture, the effect of the regeneration can be much 
more important for detectors which are sensitive 
to r's emerging at larger angles a: 30% less r's 
at 0.3 EeV for an incident neutrino flux ^{E,^) oc 
EjJ^. This is typically the case for fluorescence 
telescopes such as the ones used by HiRes (and 
Auger), which detect the faint ultra-violet light 
emitted by nitrogen molecules that are excited as 
the shower traverses the atmosphere. 

Going now to non-standard values of /3 and 
a"'^ , the combination Plow ® f^high h^deed the 
most sensitive to regeneration processes. In that 
case, neglecting the regeneration effects leads to a 
significant underestimation of the flux of emerg- 
ing taus even for incident neutrinos with rela- 
tively low energies: about 70% of the r's are lost 
when injecting neutrinos with energy 0.3 EeV. 
This statement remains true even for incident 
neutrino fluxes $(i?,/) oc E~^ when integrating 
over a small angular range (0° < a < 3°), as can 
be seen in Fig. [6l 




Figure 6. Flux of emerging taus corresponding to in- 
cident neutrinos with angle up to 3° and <I>(i5^) oc E^J^ 
using the combination (Thigh ® Plow The regeneration 
is neglected in the dashed histogram and included in the 
solid one. 
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5. CONCLUSIONS 

In our studies ^7, 8J, we have questioned the va- 
hdity of some approximations usuaUy made when 
dealing with the propagation of UHE tau leptons 
and neutrinos through the Earth. 

We have first shown that the stochastic nature 
of the radiative processes undergone by tau lep- 
tons at ultra- high energy is indeed responsible for 
large fluctuations in the tau energy losses. At the 
same time, however, these fluctuations are not 
so large as to blur the picture with respect to 
the continuous energy loss approximation, as far 
as the calculation concerns power-law injection 
spectra in a given, continuous energy range. 

We have then studied the mechanism of regen- 
eration of the Vr flux while crossing the Earth, 
to investigate its effect on the flux of emerging 
T leptons. Assuming a detector with an energy 
threshold of 0.1 EeV and sensitive only to taus 
with an emerging angle below 3°, the effect of 
regeneration is negligible for a flux of incident 
neutrinos ^{E^) oc E:^^ with the standard val- 
ues of cross-sections and energy losses. But we 
have shown that this is not valid for other as- 
sumptions on the detector performance, neither 
for less standard values of the weak cross-section 
or tau energy losses. Moreover, the error made 
in the computation of the flux of emerging taus 
will also increase with the hardness of the flux of 
incident neutrinos since the contribution of the 
regeneration effect increases with the energy of 
the incident neutrino. 

The simplification of neglecting the regenera- 
tion is thus only safe for particular values of the 
physical properties playing a role on the propaga- 
tion and specific detectors. It may lead to a sig- 
nificant underestimation of the flux of emerging 
r's when looking at non-standard values of the 
weak cross-section or tau energy losses. There- 
fore, it should be carefully treated and accounted 
for when studying the systematics due to the 
uncertainties on those properties or while using 
the Earth-skimming technique to test for instance 
higher weak cross-sections. Similarly, one should 
carefully check the effect for the characteristics of 
the actual detector before neglecting the regener- 
ation. 
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